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Abstract 

We use effective field tlieory tecliniques to compute the potentials due to spin-spin and spin- 
orbit effects, from which the 0(8182) contribution to the motion of spinning compact binaries to 
third Post-Newtonian (PN) order follow. We use a formalism which allows us to impose the spin 
supplementarity condition (SSC) in a canonical framework to all orders in the PN expansion. We 
explicitly show the equivalence with our previous results, obtained using the Newton-Wigner SSC 



at the level of the action for spin-spin and spin-orbit potentials reported in arXiv:gr-qc/0604099 
and larXiv:0712.2032[ gr-qc] respectively. 
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INTRODUCTION 



NRGR [l| , an Effective Field Theory (EFT) approach to gravity, has emerged as a pow- 
erful tool to systematically describe the dynamics of finite size objects in General Relativity 

ej -i). It has been utilized to calculate hipher order spin corrections in the PN expansion 
,l3|, dissipative effects for non-spinning [41] as well spinning objects 5|], radiation reaction 
effects in the extremal limit and corrections to thermodynamic quantities in caged black 
holes 0, B]. In this paper we will extend the formalism for spin in NRGR originally de- 
veloped in [9]. In particular we will demonstrate how to calculate the equations of motion 
(EOM) using the Routhian formalism discussed in 10|. The leading order (LO) spin-spin 
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and spin-orbit potentials were shown in to reproduce known results 
within the Newton- Wigner (NW) and covariant SSCs. In j^, the previously uncalculated 
3PN spin-spin potential was obtained using the NW SSC at the level of the action. In 10| . 
it was argued that within this approach, the Hamiltonian method is accurate up to 4PN in 
the S1S2 sector, when curvature effects start to play a role, and the canonical structure in 
the reduced phase space, (x, "P, S), is modified. However, in order to calculate the 0(8182) 
contributions to the EOM at 3PN, the spin-spin potential in {2]] is not sufficient. This fact 



was made clear by an independent calculation in 15|] where the complete potential was com- 
puted using a more traditional methodology. Within our approach, we must also include a 
correction stemming from a subleading^ffect due to spin-orbit interactions [3]. Once this 
e..a Piece is .eluded,*...... a Qa„dQ agree. As we will see, the Routhian 

formalism provides yet another independent cross check of the new results to 3PN. 

While working within the NW SSC at the level of the action is relatively simple for LO 
effects, and 0(8182) corrections up to 4PN, calculating subleading 0{Sq), or 0(8^) effects 
can be cumbersome within this methodology. Recall that the NW SSC leads to a canonical 



structure in the reduced phase space only in a fiat spacetime background [10|. For 0{Sq) 
effects this structure is lost already at 2.5PN. It is thus desirable to have a technique where 
the SSC is not imposed until the end of the calculation thus avoiding complicated algebraic 



structures. Here we elaborate upon such an approach, presented in [lO|, and compute the 
spin-spin potential to 3PN in the covariant SSC. 

Within the NRGR formalism &pin-spin, or spin-orbit, refer to the type of diagrams 
contributing to the potential 0, 9|. Since we will postpone the SSC to the later stages 
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of the computation, our results for the potentials and EOM will be written in terms of 
the spin tensor in the local frame, e.g. 5*"^. Therefore, a spin-orbit term proportional 
to S^^ in the EOM can contribute at 0(8182) once the SSC is enforced, the spin tensor 
is reduced to a three vector and the velocity in the local frame is transformed to the 
global PN frame. Furthermore, we will show that the spin-spin potential to 3PN also 
depends on S^^, and therefore it will contribute at higher orders in the 8^82 and 8i8| sectors. 



Here we will present some of the details of the calculation of the 3PN potential, as well as 
including the EOM for the spin of the constituents. We will first calculate the LO EOM due 
to the spin-orbit potential and show that we reproduce the well known results before moving 
on to the spin-spin 3PN computation. For completeness, finally we show the equivalence 
with our previous results using the NW SSC at the level of the action of , by construct- 
ing an effective potential which agrees with the results in , and from which the EOM 



follow via the canonical methods, thereby providing a formal proof of the claims in 
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io|. 



II. REVIEW OF SPIN IN GR 
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The extension of NRGR to include spin effects was achieved in |9| by adding world-line 
degrees of freedom A;^(A), which is the boost that transforms the locally flat frame (labelled 
by small Roman letters), to the co-rotating frame labelled by capital Roman letters. The 



generalized angular velocity is given by = e^*^-^^, where = e^A" and are the 
co-rotating and locally fiat basis respectively (verbeins) and e'te\,g^^ = 7]°"^. The spin S jj^u is 
introduced as the conjugate momentum to ^^y. The form of the world-line action is then 
fixed by reparametrization invariance j^, 

S = -J2{l PKd^^ + / ls'/Q\jdX,^ , (1) 

where the sum extends over the consituents , and S^"^ = S^^e^^el- Here we have not included 
higher dimensional operators which account for finite size effects. Corrections due to finite 



size effects are reported in ITj . The Mathisson-Papapetrou (MP) equations [18|] follow 
from ([1]) 0, 9|. The spin-gravity coupling in ([1]) can be rewritten by introducing the Ricci 
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rotation coefficients, uj^ = e^^D^e"-^ , as j^, [lo| 

S spin— gravity j Sj^ab^fj, U^d\^ 

(2) 

with Sf^ = S'^'^e^e^J,, the spin in the locally flat frame (we drop the L from now on). By 
further expanding ([2]) in the weak gravity limit one obtains the Feynman rules 

a a. 

Let 

us emphasize that the SSC is imposed in the local frame. 

III. THE ROUTHIAN APPROACH FOR SPINNING BODIES IN NRGR 

A Routhian formalism^ was introduced within the covariant SSC in In what follows 
we will adopt this framework and compute the 3PN corrections to the potential. The virtue 
of the Routhian formalism is that it allows us to consistently impose, and preserve upon 
evolution, the SSC in a canonical framework, and properly account for S*^ corrections to 
the potential. The price to pay is that we will work with a spin tensor, S""^, rather than a 
three vector. However, we will show later on that an effective potential in terms of (x, v, S) 
exists, which turns out to be equivalent to our previous results in 

Since the spin is a conjugate momentum, we would like to treat the spin within a Hamil- 
tonian formalism. Whereas, for the worldline position we would like to work within the 
Lagrangian formalism. That is, we would like to Legendre transform the Lagrangian with 
respect to the wordline spin degrees of freedom only. This is done within what is called the 
"Routhian" formalism 26|]. We will work in the covariant SSC, 

PaS""' = p.S^'' = 0, (3) 



with the coordinate momentum of the particle. To dynamically maintain this conditions 
we need to impose 

^(p^sn = o, (4) 



^ A similar Routhian was originally proposed in [l9| with S^^ui, = as SSC, which is equivalent to ours at 
3PN. See appendix A. 



4 



and utilizing the MP equations (which follow from ([T])) 



^^-pu pu, - ^Kp^fsb u, [b) 



yields the momentum 



P 



Notice that p ■ u = m once the SSC is imposed. 



We introduce now the following Routhian 10 1 



+ ■■■], (7) 



where the ellipses represent curvature terms necessary to account for the mismatch between 
p and M in (I3j). These terms contribute beyond the 3PN order we work in this paper^. 
In addition there are finite size corrections to ([6]) which are not shown in (JTj) but can be 
consistently included (for details see appendix A) when going to higher orders in the FN 
expansion. 

The EOM follow from 

'nd\ = o, ^ = {^«^7^}, (8) 



6x'' J dX 
where the algebra for the phase space variables {x'^,p'^, S*"*) is given by 

{x^Pj = 5^, {x^pj = 5^, {p",p/^} = 0, (9) 

{x^x'^} = 0, {p^/} = ^i^-^,fe5'^^ (10) 

{x^, 5"^} = 0, {pa, 5"''} = cu^["5^1", {P", 5"''} = (11) 

with p^ related to the canonical momentum by "P^ = p^ + ^cj^^S'''''. It can be shown that 

the EOM are equivalent to the MP equations and that the extra term in ([7]) guarantees 



In other words, to our level of accuracy, we can consider S°'^ui, = 0. 
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the preservation of the covariant SSC. Notice that with our choice of metric convention 
(+, — , — , — ), the spin vector algebra differs from the canonical S0(3) algebra by a minus sign. 
We have compensated for this convention choice by the overall minus sign in the Routhian 
of d?]) which allows us to treat TZ as the usual Lagrangian and keep the spinless feynman 
rules untouched. Therefore, the relationship between the potential and the Routhian stays 
as before for the Lagrangian, namely 

V = -7^, (13) 
and therefore the spin EOM in terms of the potential take the form, 

^Cab 



{V, S"""}. (14) 



In practice the EOM for spin can be derived from 

as one would expect, plus corrections from the 5*°* components. We will study an example 
in detail later on. 

According to the program developed in jl|, to calculate the potential we first need to 
generate a set of Feynman rules. The potential will then follow by including the appropriate 
set of Feynman diagrams. Once we have the potential in terms of the spin, position and 
velocities of the binary constituents, we can calculate the EOM using (1141) . 



A. The effective action 



Let us elaborate upon the manipulations leading to the potential. The EFT approach is 
built to separate physics at different scales. Given that the radiation and potential modes 
have a ratio of wavelengths of order v, this allows us to cleanly separate the physics of 
radiation from that of potentials in a systematic fashion. This is not to say that potential 
modes have no effect on radiation. Indeed, the tail effect arises from the coup ling of a 
radiation graviton to a potential graviton, and the EFT reproduces known results [20|| . The 
same can be said for the so-called "memory effect" . However, if we are interested in pure 
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potentials we may completely ignore the radiation mode in the effective action^. As discussed 
in l|, for spinless objects the effective NRGR action follows from the path integral {q = 1, 2) 

exp [iSNRGR[xl]] = J VH^^exp [iS[Hf,^,xl] + iScp] , (16) 

which accounts for the vacuum to vacuum amplitude in the presence of sources, in our case 
the binary. In the expression above S = Seh + Spp, that is the Einstein-Hilbert action 
plus the wordline sources, and Sgf is a suitable gauge fixing term^ jl]]. By expanding the 
instein-Hilbert action in the weak gravity limit we can immediately read off Feynman rules 



11]. Once we compute Snrgr[xI\ the EOM follow from a minimal action principle [22|], since 
we have yet to perform the path integral over the sources, namely the wordlines. Notice 
that the kinetic term is a pure phase which factors out of the path integral. Therefore, by 



summing Feynman diagrams effectively we are calculating the potential energy 



22| 



Within this framework the inclusion of spin is straightforward. The Routhian TZ will 
replace the worldline Spp action (recall TZ = —V) and the path integral in ( fT6l) will produce 
the effective potential, —Vnrgr^ from which the EOM follow via ( 1141) . Since we are not im- 
posing the SSC until the EOM is obtained, we are always dealing with a canonical structure, 
although we pay the price of having a spin tensor S*"*, rather than a three vector. The latter 
follows once the SSC imposed at the level of the EOM. For the spin dynamics we directly 
compute the potential and no kinetic piece is necessary. In a sense the spin dynamics has 
a more direct contact with the usual interpretation of the path integral as providing the 
energy of the 'vacuum' in the presence of the sources. 

Finally note that the extra terms in the action proportional to the SSC effectively act as 
Lagrange multipliers, and enforces the conservation of the SSC upon evolution. Given that 
we are dealing with second class constraints in the SSC, the Lagrange multipliers are fixed 



by preservation of the constraints 



231 ]. Moreover, since the extra term is proportional to the 



SSC itself, we are free to use the equations of motion (i.e. perform an implicit coordinate 



The LO radiation effects were computed for spinless and spinning bodies in [l| and 21 1 respectively, and 
shown to agree with known results. 
^ The gauge chosen in [1[ corresponds to an harmonic condition up to 0{G^) corrections (see Eqs. (62)-(65) 
inH). 
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shift) to replace dp/dX using (5) in the Routhian as we do in ([7]). 



IV. AN EXAMPLE: THE SPIN ORBIT EFFECTS AT LEADING ORDER 

As a warm up let us start by computing the LO effects in the EOM due to the spin-orbit 
coupling ^. In this section we will quote the contribution from each diagram to the potential. 
The relevant Feynman rules have been relegated to appendix B. Each diagrams contribution 
is of the form —iV. 

The LO spin-orbit potential is found by the instantaneous one graviton exchange diagram, 
with the LO spin vertex on one world-line and a mass vertex on the other, as discussed in 

V-,^ = (5f + Stiv'; - 2.,^)) + 1^2, (17) 

with = (g^, hence 7^?V = -V,%n- 
Applying the algebra ( fT2l) we find 

dSiii Gm2 

= [{'"'iVij - '^niV2j){6iiSkj + 6kjSii - 5ijSki - SkiSij) + {niSko - rikSm)] . (18) 
In terms of the spin 'vector', 5' = ^e'-^'^S^^, we have in the covariant SSC, 



dSi 77120 N , „ m2G 



2'-^{nx^r)xS,+ '-^{S,xn)x^r, (19) 



with V = Vi — V2. The expression in f fT9|) agrees with the known spin precession 12, [iSj 



with yU the reduced mass, after the transformation [o], [l^ 

Si = (l-^v?)Si + ivi(vi-Si). (21) 

In the expression of ( l2Ti) vi is the velocity in the local frame, which agrees with the 
coordinate veolocity in the PN (global) frame at LO. Let us add a few comments about this 



Recall that the LO spin-spin potential does not include the troublesome term proportional to 8°"^ Q . 
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distinction. First of all, we are dealing with the local spin, therefore to the order we are 
working at (recall ~ mu°- + ...), our SSC reads S°'^u\, = 0, with u"" = e^^u^. If we choose 
A = t, we have = (l,v) and what we denote as {v^,\). For the spin-spin dynamics, 
the relevant (spin-dependent) part of this relationship is (for instance for particle one) 

{;r° = 1 + (22) 

= < + ^Si'n'' + . . . , (23) 

therefore 

5f = S'^vi + 5i^e^o(xi) + . . . = (vi X Si)^ + ^ [(n x S2) x Si]^ + . . . , (24) 

where we used eg(xi) = ^(n x S2)-', which follows from the one point function, (i/g)/2, 
or simply inspection of the Kerr metric in harmonic coordinates. This will add an extra 
piece in the spin-spin EOM from the LO spin-orbit term of (ITTl) (see ( 1671) ). since there is 
modification in the algebra given by 

{Si, } = e'^'S'," = Slvi - S{v\ + ^((n X S^y S\ - (n x 82)'^^), (25) 

which will lead to a 3PN contribution in the potential^. The expression in ( l25l) is the main 
reason why we need to keep track of spin-orbit terms in the potential which will wind up 
contributing at 0(8182) in the EOM. 



A. The Equivalence of Methodologies 

1. The PN frame versus the local frame in the covariant SSC 

Naive comparison of the EOM in flTI?]) with the spin EOM in the covariant SSC, for 



instance in 



24| . shows that they are indeed different expressions. To understand the dis- 



crepancy, we have to transform from the locally flat frame where (fT9l) is defined, to the 



The extra term in (p4|) also provides an extra piece in the potential within NW SSC 
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commonly used PN frame by rotating the spin tensor using the vierbein and the metric at 
IPN order [Sli, 



s\ 



(26) 



with S''^ the spin tensor in the PN frame within the covariant SSC. One can now trace the 
disagreement back to the definition of the spin vector. In our calculations we introduced 
gjk _ ^jkiQi local frame, however, more generally we may define the spin four vector 

as 



1 



Using now fl27|) (for instance for Si) we have in the PN frame 



SI' = e 



(27) 



(28) 



Leaving 



Si = 1 + 2 



GN'm2 



1 + -^ 



r J 



Si - vi(Si ■ vi' 



(29) 



which we can expand at IPN order. 



GNm2 



Si - vi(Si ■ vi) + 



(30) 



The EOM in terms of S reads (at 1.5PN), 



dSi m2GN 
dt 



[Si(n • v) - 2n(Si • v) + (Si ■ n)(vi - 2v2)] 



which agrees with the results in 2J, |25| 



(31) 



2. Imposing the SSC before or after calculating the EOM 

It is also instructive to see how applying the SSC prior to finding the EOM leads to 
the same results as in the Routhian approach, where the SSC is enforced after the EOM 
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is obtained. In |9| it was shown how the spin EOM in the covariant SSC follow from fll7p 
once the SSC is imposed and the non canonical algebra taken into account. Let us see how 
this works for the EOM arising from the the spin-orbit interaction. Recall the commutators 
after imposing the covariant SSC are {Db stands for Dirac bracket) 



K4]o, = (32) 
= (33) 

IILq 

nab Qcd] (nO-c (J? \ nbd , f^bd Q Q \ nac 



_^^ad _ _<UL^gbc _ ^^bc _ _^^gad (34) 



with q = 1,2. (!32|) contributes a non-canoncial piece to the acceleration 



xi, 



J Db 



In this equation we have left off the other contributions in the RHS. 

Within the Routhian approach we are advocating here it is simple to show how ([35 
arises. First of all we re-write (fTTl) as 



Vr,,n = {SM + Si\v', - 2.^-)) + ^n^iSr - S^) + 1-2. (36) 

The second piece would have vanished had we imposed the covariant SSC and is thus re- 
sponsible for the new contribution on the RHS of fl35|) . It is clear that the only term which 
does not cancel out once the covariant SSC is imposed, is the one coming from -ji—^^- 
That contributes precisely the extra term in (l35l) . The resulting EOM reads (at 1.5PN) 



so _ Gn I rn2 
^1 ~ 73 



I — [-3v X Si + 6n(v X Si) ■ n + 3n ■ v(n x Si)] 



-4v X S2 + 6n(v X S2) ■ n + 6n ■ v(n x S2)} . (37) 

To establish the full equivalence we need to show that the spin equations also match. If 
we impose the covariant SSC in (!36|) and use Hamilton's equations, including the correction 
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due to fl33p . we can show that the new piece due to the algebra is (for particle one) 

SS, = ^^^(Si X n) X vi, (38) 

and the equivalence is thus proven. For more on the consistency of the Routhian approach 
see the Appendix. 



Another important point in the connection between the NW and covariant SSCs is that 
in addition to (12T|) it also entails a coordinate transformation given by 
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x,-.x,-^(v, xS,) + ..., (39) 
for g = 1, 2. Equivalently we have 

r^r-^(vxf), (40) 

2M^ sy, V J 

with ^ = + ^§2. This transformation, implemented in the LO EOM, allows us to 
transform the acceleration from the covariant ( [37] ) to the NW SSC 



V. THE 3PN SPIN SPIN POTENTIAL 

Let us now consider the 3PN potential. The result in was presented with the NW SSC 
imposed at the level of the action. Here we will derive in more detail the full expression for 
the potential in terms of the spin tensor, before imposing the covariant SSC. The resulting 
potential reproduces that of [2|] once the NW SSC is enforced. However, we will retain 
the expression in terms of S"'' and obtain the EOM via dHl). Only then will we impose 
the covariant SSC. As we will show later on, the EOM obtained with either procedure are 
equivalent in the S1S2 sector up to 4PN order as originally argued in j^. lo|. 

Following the usual rules we draw all possible Feynman diagrams which scale as v^. Each 
one of this diagrams will contribute to the effective potential by the rule —i J Vdt = diagram 
|l|, where only connected diagrams contribute. For simplicity in what follows we will 
suppress the factors of / dt. 
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To calculate the one graviton exchange contribution we should in principle draw all di- 
agrams with propagator which connects to vertices which have subleading scalings. These 
diagrams are collected in Figs. [2](a,b,c). This would be formally the correct way to do 
the calculation in the spirit of effective field theory, in that each diagram would scale ho- 
mogeneously. However, in practice it is sometimes simpler to calculate the full covariant 
one graviton diagrams and then break it into its individual pieces which scale homoge- 
neously. This allows us to calculate multiple diagrams simultaneously. The instantaneous 
one-graviton exchange can be combined into a single calculation stemming from the linear 
spin-gravity coupling. If we denote this contribution to the effective potential by Vig, and 
include also the LO piece, the combination of diagrams then reads 

-^^15= 5rV52X(^a..,/3(^l)^<xA,p(x2)) (41) 

where ( ) represents the Wick contraction. 

Let us start by considering the LO contribution from the instantaneous propagator which 
comes from the spatial components c = i,d = j, since temporal derivatives are down by v 
and Soi is down by a factor of v. The result reads 

- ^V,f = (^^) ' SfS^,^Ho.,,{x,)H,,,{x,)) (42) 

with 

dt'd^'- ^ = d,,- = \ (5,, - 3n,n,) . (43) 

|xi - X2I J r ^ ^ 

From here we read off the LO spin-spin potential (see Fig. [1]) 



« = -^(si-S,-3i:^), (44) 



Now let us consider the subleading contributions. There are multiple terms at 3PN. A 
factor of V arises from either a spatio-temporal component S'oi, a temporal derivative, or an 
explicit factor of v. Expanding out and keeping only the terms which contribute at 3PN 
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T 



FIG. 1: Leading order spin-spin interaction. 



gives 



2M, 



pi 



Sf {HoifiHoi^k) + 'S'i*5'™"'yf (iffco,j-f^om,n) + 1^2)] 



(45) 




V 



-o- 



V 

c) 



FIG. 2: Diagrams contributing to 3PN order which do not involve non-linear interactions. The 
blob represents a spin insertion and the cross corresponds to a propagator correction. 



V 

-O 



V 

-o 



n 2 

V V 

a) 



y2 ^0 



FIG. 3: Non-linear contributions to the 3PN spin-spin potential. 



The evaluation of these integrals is straightforward. Here we evaluate one particular 
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integal which needs more dehcate consideration. Namely the contribution where one v 
comes from a temporal derivative while the other comes from 5**^* . This contribution is given 
by 

= _tsfS^2 [ ^"^^^^ (-g2^-»poai(Ai)-t2(A2))^^g2^-»p-(xi(Ai)-X2(A2))^ 

2 J P 

(46) 

Now recall that_^ ~ v'^^S^ so we can neglect the time variation of spin at the order we 
are working at |9|. It is therefore convenient to trade dt^ for dt^ picking up a minus sign, 
and integrating by parts with no net effect. Had we kept dt^ we would have to deal with 
^5"°* which we can not neglect. To see this notice that, imposing the SSC would introduce 
an acceleration dependent piece into the Lagrangian, which can be eliminated using the LO 



EOM. However doing so entails a change of coordinates 



291] which is not preferable^ . If we 



consider now the instantaneous interaction once again we get 

S'^Si'{Ho.,o{x^)Hoi,k{x2)) = 5(ti -t2)^^f^f5fc'5o^y"i^(e-^P-(^^('^)-^^('^))) 

= m - t2)^5f^f (Sn'^n ■ V2 - v^^)- (47) 

The net result from the instantaneous one graviton exchange is then given by 

- iVi;;'' = '-^{6'^ - 3nV) [^f + vi ■ ^r2S^St + ^^^'2'^^ '^2" (48) 
- y'lv^Si'^Sr + Sf'Stiv!^ - v'^) + 5f (t;f - vl^) 

+ (3n'v2 ■ n - vl)S°,^Sl' + (3n'vi ■ n - v[)S^^S'['] (49) 

The first corrections to instantaneity comes from the diagram shown in Fig. [2t. This 
correction comes from expanding the graviton propapagtor to second order in the energy 



^ Incidentally, had we insisted on keeping dt^ and included this acceleration piece, it turns out that imposing 
the NW SSC at the level of the action would reproduce the exact same form for the 3PN Hamiltonian 
of 



15|. That is actually the case due to a cancelation of this 0{G^) acceleration piece with the extra 
term stemming from subleading corrections in the spin-orbit potential due to (j25p . A similar result can 
be foimd in 30]. Within the Routhian approach the trade for dt2 is preferable. 
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which is suppressed by a factor of v relative to the spatial momentum, i.e 

Po - P^ P^ 

The result from this diagram is very similar to the spinless case and we have 

FigEfc = -i^S^'S^' [vi ■ V2(5*^' - Sn'n^) - 3vi ■ nvs ■ n(5'' - 5n'n^) 

— 3v2 ■ n{v{n' + vln^) — 3vi ■ n{vln'' + v\n^) + v\v{ + v\v{\ (51) 

Let us now consider the terms involving non-linear graviton interactions, as shown in 
Figs. [3fc. and[3]D. Let us start with the so called "seagull" topology 

Fig Wp = ^S'rs'i {T{Hoi,^{x^)H^,{x{){HjHox,i + Ho,,kHkj))). (52) 

Note that for this diagram there are a few Wick contractions, that is, the two graviton 
vertex can contract in two ways with the mass and spin vertex on the opposing world line. 
However many of these contractions vanish since index structures vanish, e.g. 

{HooHoi) = 0. (53) 

The result is given by 

FigEb = ^^^^i^Sf ^^^(SnV - 26''') (54) 

Finally we have the diagram with the three graviton interaction. Again there are mutiple 
ways of doing the Wick contractions. As previously mentioned, the best way to handle 



these contraction is using a symbolic manipulation program j27|] where symmetrization is 
simply handled. The integrals are all variations of the same result used in the one graviton 
exchange diagram. The result for the diagram is 



Fig m = tmi^Stsi'{4:n''n' - 



(55) 
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Note that in calculating this diagram we encounter multiple power divergent integrals. 
These divergences can be absorbed into pure counter-terms, since they just renormalize 
the mass and possibly other quadrupole moments. It is simple to see how this occurs in a 
diagramatic language. The divergences occur when one of the propagators ending on the 
line which has a mass insertion is cancelled by a power of arises from the momentum 
depedence of the three graviton vertex. One of the lines in the diagrams then contracts 
to a point. The resulting diagram looks like an interaction between a self energy (mass 
correction) and the spin on the opposite line. The result of this renormalization is that we 
may simply drop these divergent integrals. As was explained in jl| no physical logarithmic 
divergences occur until 5PN order for the case of spinless particles. As it was shown in 9| 
that is also the case for spinning bodies and logarithmic divergences due to finite size effects 
do no show up until 0{v^^). This generalizes the so called "effacement" of internal structure 



32l | to the case of spinning bodies j^. The logarithmic divergences are renormalized by 
absorption into finite size parameters which present a non trivial renormalization group 
flow. These are tidally induced effects which in turn do not contribute to the metric 
solution as it is expected from Birkhoff 's theorem. However, there are other types of finite 
size effects, the so called self-induced effects, which do appear at lower orders as explained 
in 91]. This kind of effects are encoded in operators whose coefficients are fixed, like the 
mass, and can be generated by power law divergences 9|]. For instance in the case of a 
rotating black hole, finite size corrections appear due to the quadrupole moment of the Kerr 



spacetime. The coefficient is set by the Kerr metric and it is proportional to 5"^ 3l|]. The 
LO corrections (at 2PN) were computed in j^, and subleading effects are reported in ~v\ . 

Gathering all the pieces together, plus mirror images, we have the complete spin-spin 
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potential to 3PN, prior to imposing the covariant SSC 
Gn 



V 



spin 



+ SfSi\v^^ - vl) + S\^Si\v\ - t;2^)) + ]^Sl'S''J (3vi ■ nva ■ - ^n'n^) 
+ 3vi ■ n{v{n'' + v\n^) + 3v2 ■ n(wjn* + v\n^) - v\v{ - t^i^j) 
+ (3n'v2 ■ n - v[)SfSf + (3n'vi ■ n - v[)SfSl^] 
+ {^- S{^Si\5^^ - ?>n^n^) 

+ ^n^- (5f + - H)) - (5f + Si\vl - 2v\)) , (56) 



where we included the LO spin-orbit term which will be relevant latter on due to flMj) and 
( l25l) . The spin potential in (l56l) is the main result of the paper from which the EOM to 3PN 
order can obtained via ffUl) . 



A. The spin Hamiltonian in the NW SSC to 3PN 

Notice that the spin-spin part of the expression in fl56l) agrees with the result reported 
in ^ if we impose the NW SSC^. However, as we mentioned earlier, to obtain all the 
contribution in the S1S2 sector we need to include subleading corrections in the spin-orbit 
potential coming from (l24l) . The extra term takes the form 3| 

^ (m2n'5re^o(xi) - rmn^Sl^eii^,)) = ^ ((Si x n) ■ (n x S2)) . (57) 



It also provides the spin-orbit potential from which the precession equation follows 
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For completeness, we present the S1S2 potential in the NW SSC to 3PN order 



V, 



sls2 
NW 



G 



N 



2r3 



Si ■ S2 ^^vi ■ V2 - 3vi ■ nv2 ■ n - (v^ + Vg)^ - Si ■ V1S2 ■ V2 



- -Si ■ V2S2 ■ vi + Si ■ V2S2 ■ V2 + S2 ■ viSi ■ vi + 3Si ■ nS2 ■ n (vi ■ V2 + 5vi ■ nv2 ■ n) 

— 3Si ■ V1S2 ■ nv2 • n — 3S2 ■ V2S1 ■ nvi ■ n + 3(v2 x Si) ■ n(v2 x S2) ■ n 

3 

+ 3(vi X Si) • n(vi X S2) ■ n - -(vi x Si) ■ n(v2 x S2) ■ n - 6(vi x S2) ■ n(v2 x Si) ■ n 



2r4 



(5Si ■ S2 — 17Si ■ nS2 ■ n) 



G 



N 



(Si ■ S2 - 3Si ■ nS2 ■ n) 



(58) 



As it was argued in [10( the EOM in the S1S2 sector follow from flSSl) by means of the 
'traditional' Hamiltonian approach up to 4PN order 3|]. The spin dependent part of the 
Hamiltonian can be readily obtained from (158|) . (14^ and (fT71) . and takes the form (ignoring 
2.5PN spin-orbit and 3PN spin^ terms) to 3PN 



TTSpm 



G 



N 



-iVi X Si) ■ n(P2 X S2) ■ n + 6(^2 x Si) ■ n(Pi x S2) ■ n 



(59) 



-15(Pi ■ n)(P2 ■ n)(Si ■ n)(S2 ■ n) + ^{V^ ■ Si)(Pi ■ S2) - ^(^2 ■ Pi)(Si ■ S2) 
-3(Pi ■ P2)(Si • n)(S2 ■ n) + 3(Pi ■ Si)(n ■ P2)(S2 ■ n) + 3(^2 ■ S2)(n ■ Pi)(Si ■ n) 



-3(^2 ■ n)(n ■ Pi)(Si ■ S2) + (V2 ■ S2)(Pi ■ Si)] 

[VliS^ ■ S2) - 3(Pi X Si) ■ n(Pi X S2) ■ n - (Pi ■ S2)(Pi ■ Si)] 



[P|(Si ■ S2) - 3(^2 X Si) ■ n(P2 X S2) ■ n - {V2 ■ S2)(P2 ■ Si)] 
^ G^(mi + m2) ^^^g^ _ _ 23(g^ . ^^^g^ " n)) - ^ (Si ■ S2 - 3(Si ■ n)(S2 ■ n)) 



G 



N 



2r4 

^""^ X Vi) ■ Si - 2(n x V2) ■ Si + 2(n x Vi) ■ S2 - |^(n x V2) ■ S2 



2mi 



2m2 



where 



Vi = mivi + 2 —n x S2 H ^^5— n x Si, 1 ^ 2. 



(60) 



To obtain the EOM however we will proceed differently, and again we will not impose the 
SSC up until after we have solved for the EOM resulting from the potential in fl56|) using the 
Routhian approach. As we shall see the extra piece due to spin-orbit effects will come from 
( |25|1 . We will explicitly show however that the results are equivalent. Let us remark that a 
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Hamiltonian similar to that in (M?]) was recently found in 15|, and shown to be equivalent 
in [sl once the spin-orbit effect is included. 



VI. THE SPIN EQUATION OF MOTION TO 3PN ORDER 

The 3PN contribution to the EOM for spin follows from the potential in a similar fashion 
to the LO spin-orbit example. Let us proceed systematically for particle one. For the 
spin-spin part of the potential to 3PN in (!56|) we have two pieces, one depending on Si, 



G 



N 



Si ■ S2 ^^vi ■ V2 - V2 - ■ vin ■ + ■ 8211 ■ Si (vi ■ V2 + 5n ■ V2n ■ Vi) 

+Si ■ V2S2 ■ V2 - -^Si ■ V1S2 ■ V2 - ^S2 ■ viSi ■ V2 - 3n ■ (v2 X Si)n ■ (vi x S2) 

3 3 
+3n ■ (v2 X Si)n ■ (v2 x S2) + -n ■ Vin ■ S2S1 ■ V2 - -n ■ Vin ■ S1S2 ■ V2 

3 3 
- -n ■ V2n • S1S2 ■ vi - -n ■ V2n ■ S2S1 ■ Vi 



+ 



Gn SMGl 



+ 



N 



[Si ■ S2 - 3Si ■ nS2 ■ n) 



(61) 



and another one 



Ai 



Gn 

{(3v2 - vi) X S2 - 3n ■ (2v2 - Vi) x S2)n - 3n ■ V2(n x S2)} 



G 



a!(2) + (2vi X S2 - 3n ■ (vi x S2)n - 3n ■ vi(n x S2)) 



(62) 



(63) 



with a^2) the S2 part of the acceleration in the local frame. The latter is given by 



d fG 



(64) 



where is the acceleration in the PN frame given in (1371) . Then (with x = S2 + ^Si^ 



G 



N 



[— 3v X ;^ + 6n(v x ;^) ■ n + 3n ■ v(n x x)] ■ 



(65) 
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Notice we also have 



Ai = Sil%) + vi X w,"' 



) 



(66) 



with cUq* the LO spin-spin frequency. This expression will be useful later on to prove the 
equivalence with our previous results in 

Using ([Hj) the 0(8182) part of the spin EOM ends up being 



dSi 



{uo'q' + cuf ) X 81 + (vi X 81) X Ai 



N 



n X [(n X 82) X 8 



(67) 



where the last term follows from the correction in the spin-orbit part of the potential in 
(|56|) due to (125]), and 



G 



N 



G 



N 



82 - 3n82 ■ n) (68) 

/I 3 \ 3 

82 ( -vi ■ V2 - V2 - -n ■ vin ■ V2 j + 2^(^ ' '^2) (vi ■ V2 + 5n ■ V2n ■ Vi) 

+V2(82 • V2) - ^vi(82 ■ V2) - ^V2(82 ■ vi) - 3(n x V2)n ■ (vi x 82) 

3 3 
+3(n X V2)n ■ (v2 x 82) + -V2(n ■ Vi)(n ■ 82) - 2^i^ ' vi)(82 ■ V2) 



3 3 

-n(n ■ V2)(82 ■ vi) - (^^ ' V2)(n ■ 82) 



3MG 



N 



3n82 ■ n) 



(69) 



In what follows we will show how to reproduce the precession equation and the equivalence 
with the result of 



VII. THE PRECESSION EQUATION TO 3PN AND THE EQUIVALENCE WITH 
OUR PREVIOUS RESULTS USING THE NW SSC 



The precession equation, 



dS 



^ = ^7 X S., (70) 
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can be obtained from (l60l) with oj, 



dS„ 



[q = 1,2), for instance for particle 1 (ignoring 



linear in spin, and also spin^, terms). 



Gn 
2^ 



-n X Vi(v2 X S2) • n + 6n X V2(vi x S2) ■ n 



3 3 
-15n(vi ■ n)(v2 ■ n)(S2 ■ n) + -V2(vi ■ S2) - -(v2 ■ Vi)S2 

-3n(vi ■ V2)(S2 ■ n) + 3vi(n ■ V2)(S2 ■ n) + 3n(v2 ■ S2)(n ■ Vi) 

+3S2(v2 ■ n)(n ■ Vi) + Vi(v2 ■ S2)] 



2r3 
Gn 



[v^S2 - 3n X vi(vi X S2) ■ n - Vi(vi • S2)] 
[V2S2 - 3n X V2(v2 X S2) ■ n - V2(v2 ■ S2)] 



Gifmi + 7712) 



G 



(5S2 - 17n(S2 ■ n)) - ^ (S2 - 3n(S2 ■ n)) 



2r4 



(71) 



Notice that equivalently we have ujg 



C)T/Sls2 

%f(? = l>2). 



In what follows we will show (170|) is equivalent to (167|) up to O^Sf) effects. To transform 
the EOM in covariant SSC to NW SSC, as a first step we need to implement the shifts in 
( !2T|) and (!39|) . Recall (12T!) already transforms the LO spin-orbit part of the EOM into a 
precession equation (see (|T9|) and ( l20l) ). Also the coordinate transformation in ( l39i) shifts 
the form of the frequency in the precession equation from the spin-orbit part. The EOM in 
terms for Si reads 



with 



dSi 
~dt 



cDf X Si 



5ul° + u'q' + 5u'q' + cjf + -vi X Ai + - 



1 m2G 



[So X n) X n 



(72) 



(73) 



r~so Gn 



n X -vi 
' 2 



6vc 



1712 



— (n X vi) ■ Si - (n X V2) ■ S2 

nil 



m2 



+ V2 X S2 Vi X Si 

nil 



X 



-Vi 

2 



2v, 



G%m2 



n X 



3m2 
4mi 



:n X Sil 



nil I 
ni2 



n X S2 



(74) 
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and 



G 



N 



2^3 



V2S2 - V2(S2 ■ V2) - 3n(n ■ S2)v2 + 3n(n ■ V2)(v2 ■ S2) 



(75) 



From the expression in fl74p we will only consider the S1S2 contributions. What we need 
now is to find an additional contributions (curvature effects) which would take the form of 
fl73|) into the expression in (17T1) . First of all notice that Ai ends up effectively like in the NW 
SSC, due to the | in (!73|) . However, there is a piece which differs from the full NW form 
and comes from the 5f 5f term in the potential. For the expression in cuf^ the difference is 
just the factor of | for S^^ in the NW SSC. Henceforth, we can split the terms in (1731) as 



Ai 



G 



Ar - ^ ( 3n(n x V2) ■ S2 - V2 x S2 



.w _ ^ 

+3n X (v2 X S2)(n ■ vi 



V2 X S2) X (2vi - V2) - 3n X (vi - V2)(n x V2) ■ S2 



(76) 



(77) 



Notice that u'^'^ = ^vi x A^"" + u'^'" and the EOM becomes 



dSi dSi 
~~dt ^ 



dt 



Gn 
2^ 



-3(n X V2)(n x V2) ■ S2 + (v2 x S2) x (vi + V2) + V2S2 



+3n X (v2 X S2)(n ■ Vi) - V2(S2 ■ V2) - 3n(n ■ S2)v2 + 3n(n ■ V2)(v2 ■ S2) 
G%mi 



2r^ 



n X n X S-: 



X Si 



X Si 
(78) 



The extra shift we need to add to fl21l) and transform away the undesired pieces ends up 
being 



Sr = Si + ||(S2(v2-n) 



(S2 ■ n)v2 X S 



G 



;i - ^v?)Si + -vi(v ■ Si) + [(V2 X S2) X n] X Si + . 



(79) 
(80) 



and the equivalence is thus formally proven. Notice that S""' = S + 0{Gn). Recall that 
S reproduces the spin dynamics in NW gauge at LO, in particular the LO precession 
equation in (l20l) . However, at next to LO, to transform to the NW gauge we needed to 
take into account curvature effects that modify the shift in ( l2Ti) . Some of these effects 
were already included once the velocity in the local frame is transformed to the coordinate 
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velocity in the PN frame (see fl2ip and (1231] ). The extra term, necessary to account for 
the discrepancy between a flat and curved background, appears in flHOl) . Notice that in 
the limit Gjv — these contributions vanish. To avoid confusion, to 3PN one can skip 
the intermediate step in fl2T|) which defines S, and use flHOj) to relate S with S""", the 
local spin in the covariant and NW SSC respectively. The equivalence of results thus follows. 

To show that the position dynamics is also recovered, once the spin EOM is reproduced, 
we can simply construct an effective potential as 

Veff = (sr) ■ sr (81) 

from which the spin corrections to the position dynamics can be derived via the 'traditional 



Hamiltonian approach. The above expression obviously reproduces the results of 
Nevertheless, the more traditional spin dynamics in covariant SSC is shown in (167|) with the 
spin defined in the local frame. To transform to the PN frame one can proceed as we did in 
([26])- (EI]) for the LO case. 

A. Adding terms 

As we mentioned earlier, the spin EOM which follows from ( l58l) fails to reproduce all of 
the 0{Sg) terms. These terms can be computed by working within the Routhian formalism 
by adding the corrections due to the Riemann dependent term in We may use this term 
as written or equivalently we may perform a field redefinition such that 

ndeab^ ^ ^= — TT !==-■ (S^j 



The procedure for calculating the potential follows the exact same steps as before and it 
can be shown that 



17j the potential due to this term in the covariant SSC takes the form 
F/;^ = ...-(a?^,))'5S, + 1^2, (83) 
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from which we get the following contribution to the spin EOM 



dSi 
~dt 



...+ (al^i) X Si) X vi 



(84) 



with a^j^^ the Si dependent part of the spin-orbit acceleration in the local frame (see ( l65l) ). 
This is however not yet complete since we are still missing 5*^ corrections stemming from 
finite size effects, as well as non-linear corrections (~ G\) as in Fig. [3K, with two LO spin 
insertions on the same worldline^. Finite size effects on the other hand are encoded in higher 
dimensional operators jl, 9]. For the case of self-induced spin effects the new term in the 
Routhian takes the form (g = 1, 2) 



^£52 a ^cb 



(85) 



in the worldline 



9], where S'^'' is defined as[l 



^ab ^ gab ^ !f£5c[a^6] 



(86) 



which guarantees the SSC is preserved in time (see appendix). In the expression of fl85l) the 
Wilson coefficients, Cl^L, are constants which are determined solely by the nature of the 

riQ 

object 19|, and Eab is the electric component of the Weyl tensor in the local frame. In the 
case of a rotating black hole we have Ces^ = ^"^^ this term represents the non- vanishing 
quadrupole moment of the Kerr solution. The LO self-induced finite size contribution to 
the potential thus takes the form 

m 



2PN 



-c 



(1) rn2 



ES^ 2mir3 



fSi ■ Si - 3Si ■ nSi ■ n) + 1 



(87) 



Higher order corrections will follow from fl85]) in a similar manner. We report the full 
0{Sg) contribution in a companion publication ITj. 



^ Notice that the would be 3PN contribution from a seagull diagram similar to Fig. [Sh, where a non-linear 
contribution from the mass worldline couples to two spin insertions in the companion worldline. vanishes. 
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VIII. CONCLUSIONS 



In this paper we have presented the details of the calculation of the 0(8182) effects to 
3PN order. We computed the potential, and showed how to calculate using a Routhian 
approach, imposing the SSC only at the last stage of the calculation. The EOM follow from 
fl56|) via f|T4|) . We proved the equivalence of this methodology, with a covariant SSC, to 



that originally espoused in 



Sj, where we calculated within the NW imposing the SSC at 
the level of the action. In this paper we have not included effects which go as S"^ such as 
:inite size effects. The first non-zero finite size effects for spinless particles start out at 5PN 







m, whereas spin induced fin.te si.e effects show up at LO, e^g^ 2PN «. TidaUy induced 

finite size effects (logarithmic effects) first appear at 5PN for the case of spinning bodies [9| . 
In a subsequent paper we will present the next to LO 5^ effects using the same formalism 
discussed here. 

This work was supported in part by the Department of Energy under Grants DOE-ER- 
40682-143 and DEAC02-6CH03000. RAP also acknowledges support from the Foundational 
Questions Institute (fqxi.org) under grant RPFI-06-18, and funds from the University of 
California. 
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APPENDIX A: ON THE ROUTHIAN FORMALISM 



Here we discuss some subtleties of the Routhian formalism, in particular its consistency 
with regards of the preservation upon time evolution of the SSC and the equivalence 
between imposing the SSC before or after obtaining the EOM. 



Let us start with the expression in (l82l) . This is nothing but performing a field redefinition 
|l| (or coordinate transformation) given by 5x''(A) ~ S'^'^Uc in the worldline action. Notice 
that it vanishes when one imposes the SSC. This implies that we could have indeed started 
using this other form if we wished, since the MP equations are also recovered. Notice that 
this extra acceleration dependent piece effectively entails adding a term^° in the Routhian, 

^^V^, (Al) 

and a modified gravity-spin interaction of the form [1 

- (A2) 



with 



gab ^ gab j^^gcia^h]^ (A3) 



V? 



The term in flAip becomes crucial, and generates a piece iidS'^Uc into the potential (recall 
TZ = —V). It is simple to show that these extra terms do not effect the LO spin potentials, 
and that the 3PN results reported in this paper are also reproduced. Notice we have now 
S"''^Uf, = algebraically. Written this way the SSC is manifestly preserved as we will now 
show. The algebra in terms of is that of ([12]) with t]"^ r/"^ - ^ as in ([34]). We can 
then show that 

Ub{S''\S''^} = 0. (A4) 

From here we have, using ([8]), 



-XS'^^'ub) = Ub{S''\ TZoiS'^')} + Ud^{S''\ S'^'jut + S'^'iib = 0, (A5) 



Recall from (O that p'^ ^ m-^. 



27 



since 

«45'^^7^o(5"')} = 0, (A6) 

where TZq stands for the Routhian without the acceleration dependent term. The same 
obviously follows from the Routhian in (JTj). The expression in ( ]A6[) guarantees that higher 
dimensional operators written in terms of S"'^ will preserve the SSC upon evolution. 



Regarding the equivalence between imposing the SSC before or after obtaining the EOM, 
we will show that the extra piece in the equations of motion in the Routhian formulation 
exactly reproduces the Dirac bracket structure when imposing the SSC at the level of the 
action. We start by noticing that the LO Dirac algebra of Sab agrees with that of Poisson 
algebra of iS"*. Since the only term in the Routhian which does not depend on S"'^ is the 
one in flAip . the above equivalence seems to rely on whether flAll) can account for the extra 
pieces induced by the non-canonical algebra of fl32|33|34l) . In what follows we show that is 
the case. The extra acceleration piece in the Routhian produces a term in the spin EOM 
given by (at LO (a° ~ 0) ) 

{^"^ S"^}udUc ~ Uc {S'"'u'' - S'^'u") . (A7) 
After imposing the covariant SSC (using Dirac brackets), ^{Sab'uf') = leads to 

Uh[S''',H{x,p)]Db = S'''uh, (A8) 



which implies 

J cab 

= H{x, p)]Db = UciS^'u" - S'^u-) + ( A9) 

with F"'' = -F^" and F^^ui, = 0. By comparison with flA6P we immediately recognize that 
F""^ is nothing but the Poisson bracket structure in [ , Therefore the extra piece in 

(IA9I) comes from the non-canonical part of the Dirac bracket as we advertised above, which 
agrees with the acceleration dependent term in (]A7p . and the equivalence follows. 
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k 




m/mp 


1/r 







TABLE I: NRGR power counting rules for potential modes. 



APPENDIX B: FEYNMAN RULES: SPIN GRAVITON VERTEX 

The spin-less part of the Feynman rules are identical to those in The potential 
graviton propagator in the gauge of [l| is given by 



d^k 1 



ikx 



(Bl) 



where 



The mass vertices can be read off from TZ 



(B2) 



a=l,2 P 



2-^00 + HoiVai + -Hoq\^ + -HijVaiVaj 



+ ... 



(B3) 



where we have only included terms which are suppressed by v'^, as higher order terms will 
not contribute at 3PN once full diagrams are computed. In this expression, and from now 
on, we chose X = = t, and therefore = (1, ^ = v). The fields have arguments which 
are the wordline coordinates and there is an implied affine parameter integral. Each dia- 
gram will contain an overall time integration which is dropped when extracting the potential. 



Following standard power counting procedures one arrives at the scaling laws for the 
NRGR fields shown in table I. In the last column we have introduced m'i = V. the Planck 
mass and L = mvr the angular momentum, with v, r the relative velocity and orbit scale 
respectively. We will not consider radiation in this paper. By including the appropriate spin 
vertices the higher order radiation can be calculaled following the methodology introduced 
in Q. 

The non-linear graviton interactions are obtained from expanding out the Einstein- 
Hilbert action The Feynman rule for the three graviton vertex has not been written 
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down due to its length^^. A Mathematica code showing how to include this vertex can be 
found at (27|. 



In order to obtain the spin-graviton vertex we also need to expand the metric in the weak 
gravity limit. In terms of the verbein we have 



p"-p'> n u = n + = 4- Sp" Sp" = — !— ff" —H^W + fB4l 



Using (lB4p we can now expand the __Ricci coefficients u'^^ in the weak gravity limit and 
extract the spin-graviton vertex rules 



tNRGR _ _J_TT qik f-nf-\ 

^™ = ^ {H^,,kS"'v^ + Hoo,kS'') , (B6) 

tNRGR _ I TT COk„J , TT QiO\ 

^ S'^{H}Hox,^-H!;Ho,,k). (B7) 



4m2 



The appropriated scaling of each vertex is derived from the rules in ([T]). Furthermore, 
the components of the spin tensor scale as 5*°* ~ v^S^\ and S ^ Lv for maximally rotating 
compact bodies. 
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